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Abstract 

The cross section for e + e~ — > 7r + 7r~^/>(2S) between threshold and y/s = 5.5 GeV is measured 
using 673 fb^ 1 of data on and off the T(4S) resonance collected with the Belle detector at KEKB. 
Two resonant structures are observed in the Tr + Tr~ifj(2S) invariant mass distribution, one at 4361 ± 
9 ± 9 MeV/c 2 with a width of 74 ± 15 ± 10 MeV/c 2 , and another at 4664 ± 11 ± 5 MeV/c 2 with a 
width of 48 ± 15 ± 3 MeV/c 2 , if the mass spectrum is parameterized with the coherent sum of two 
Breit-Wigner functions. These values do not match those of any of the known charmonium states. 

PACS numbers: 14.40. Gx, 13.25.Gv, 13.66.Bc 
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In a recently reported study of the initial state radiation (ISR) process, e + e" — > 
7/5H 7r+7r ~J/ipi the BaBar Collaboration observed an accumulation of events near 
4.26 GeV/c 2 in the invariant-mass spectrum of 7r + 7T~J/ip [l[ that they attributed to a 
possible new resonance, the F(4260). This structure was also observed by the CLEO ^ 
and Belle Collaborations using the same technique 0; in addition, there is a broad struc- 
ture near 4.05 GeV/c 2 in the Belle data. In a subsequent search for the y(4260) in the 
e + e~ — > 7/5R7r + 7r^^(25') process, the BaBar Collaboration observed a different structure 
at m = 4324 ± 24 MeV/c 2 with a width of 172 ± 33 MeV/c 2 fl that is neither consistent 
with the y(4260) — > 7T + 7r~ip(2S) peak nor with ■0(4415) — > n + Tr~ip(2S) decay. There are 
now more observed J PC = 1 states than predicted by potential models 0] in the mass 
region between 3.8 GeV/c 2 and 4.5 GeV/c 2 ; it is possible that one or more of these new 
states are exotic. However, it should be noted that other interpretations that do not require 
resonances have been proposed 0. 



In this Letter, we report an investigation of the e + e _ — > n + ii~il)(2S) process using ISR 
events observed with the Belle detector [?J at the KEKB asymmetric-energy e + e~ (3.5 on 
8 GeV) collider Here ip{2S) is reconstructed in the 7t + tt~ J/ijj — > Tc + 7r~£ + £~ (I = e,/x) 
final state. The integrated luminosity used in this analysis is 673 fb" 1 . About 90% of the 
data were collected at the T(AS) resonance (y/s = 10.58 GeV), and the rest were taken at 
a center-of-mass (CM) energy 60 MeV below the T(45) peak. 

We use the PHOKHARA event generator [9J to simulate the process e + e~ — > 
1 1 sr^ k~ %p{2S) . In the generator, one or two photons may be emitted before forming the 
resonance X, which then decays to 7T + 7r~ip(2S), with ip(2S) — > n + ix~ J / ip and J/ip — > e + e~ 
or When generating X —> ii + ix~ip{2S), a pure S'-wave between the ttti system and the 

i/j(2S), as well as between the 7r + and 7r~ is assumed. The kinematics of X decays are mod- 
elled with the 7T7T invariant mass distribution observed in our data, while ip(2S) — > n + it~ J / ip 
events are generated according to previous measurements T3 • 



For a candidate event, we require six good charged tracks with zero net charge. A good 
charged track has transverse momentum greater than 0.1 GeV/c and impact parameters 
with respect to the interaction point of dr < 0.5 cm in the r-(j) plane and \dz\ < 5 (2) cm in 
the r-z plane for pions (leptons). For each charged track, information from different detector 



subsystems is combined to form a likelihood for each particle species (i), Tracks 



with TZk = r C fr < 0-4, are identified as pions with an efficiency of about 95% for the 
tracks of interest. Similar likelihood ratios are formed for electron and muon identification. 
For electrons from J/ip e + e~, both tracks are required to have lZ e > 0.1. For muons 
from J/ijj — > fi + fi~, one of the tracks is required to have TZ^ > 0.95; in addition, if one 
of the muon candidates has no muon identification (ID) information, the polar angles of 
the two muon candidates in the 7r + 7r~/i + /i~ center-of-mass system are required to satisfy 
I cos# At | < 0.75, based on a comparison between data and MC simulation. The lepton ID 
efficiency is about 90% for J/ip — > e + e~ and 87% for J/ijj — > fi + ^r. The detection of the 
ISR photon is not required; instead, we require |M 2 CC | < 2.0 (GeV/c 2 ) 2 , where M 2 ec is the 
square of the mass recoiling against the six charged particle system assuming that four of 
them are pions and the other two are either electrons or muons. Events with 7-conversions 
are removed by requiring TZ e < 0.75 for the tt + tt~ tracks accompanying the if)(2S). 

The dilepton invariant mass distribution (the bremsstrahlung photons in the e + e~ fi- 
nal state are included) for events that survive these selection requirements is shown in 
Fig. [T](a); it is fitted with a Gaussian and a second-order polynomial. A dilepton pair is 
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considered as a J/ip candidate if its invariant mass (me+i-) is within ±45 MeV/c 2 (the 
mass resolution is 16 MeV/c 2 ) of the J/ip nominal mass (mj/^). If there are multiple 
7r + 7r~ combinations that satisfy the ip{2S) requirements, the one with \m-n+-n-i+i- — mg+£-\, 
the mass difference between the ip{2S) and J/ip [1^ . closest to 0.589 GeV/c 2 is selected; 
here m 7r + 7r -^+.£- is the invariant mass of the 7r + 7T~£ + £~ system. Figure QJb) shows the 
m^+^-j/^ (= m 7r + 7r -£+^- — mi+i- + mj/^) distribution. Fitting with a Gaussian and a 
second-order polynomial yields a mass resolution of 3 MeV/c 2 . We define a ?p(2S) sig- 
nal region as m^+^-j/^ G [3.67,3.70] GeV/c 2 , and a ip(2S) mass sideband region as 
m n + n -j/^ G [3.64,3.67] GeV/c 2 or m^+^-j/^ G [3.70,3.73] GeV/c 2 , which is twice as wide 
as the signal region. 




2.9 3 3.1 3.2 3.3 3.64 3.69 3.74 
M(l + r) (GeV/c 2 ) M(7cVJA|/) (GeV/c 2 ) 

FIG. 1: Invariant mass distributions of £ + £~ (a) and ir + ir~ J/tp (b) for selected tt + tt~ tt + tt~ £ + £~ 
candidates. The curves show fits described in the text. 

Figure [2] shows the 7i + n~ip(2S) invariant mass ("7,^+^-^(25) = mn+n-Tr+n-e+e- — 
m^+^i+i- + 771^,(25), where m^s) is the nominal ?p(2S) mass) for selected i/j(2S) events, 
together with background estimated from the scaled ip(2S) mass sidebands. Two distinct 
peaks are evident in Fig. [21 one at 4.36 GeV/c 2 and another at 4.66 GeV/c 2 . As can be 
seen from the plot, the background determined from the ip(2S) mass sidebands is very low. 
Backgrounds not described by the sidebands are negligible; these include n + 7i~ip(2S) events, 
in which the ip(2S) does not decay to n+n^J/ip (J/ip — > £ + £~), and events with a ip{2S) 
and other particles instead of 7r + n~ in the final state. 

Figure [3] shows the M 2 ec and polar angle distributions of the ir + 7r~ip(2S) system in the 
e + e~ CM frame for n + 'n~ip{2S) events with m 7r + 7r -^( 2 5) G [4.0,5.5] GeV/c 2 . The data agree 
with the MC simulation (shown as histograms) well, indicating that the signal events are 
produced via I SR. Figure H] shows the 7t + tt~ invariant mass distributions for events with 
m 7r + 7r -^(25) G [4.0,4.5] GeV/c 2 , and m n + n -^^ 2 s) G [4.5,4.9] GeV/c 2 . In both cases, the mass 
distributions differ from the phase-space expectation and tend to be concentrated at high 
mass. In the high mass resonance region, most of the tt + tt~ candidates are consistent with 
a / (980) decay. 

An unbinned maximum likelihood fit that includes two coherent P-wave Breit-Wigner 
(BW) functions and a constant, incoherent background is applied to the 7r + 7r _ ^(2S') mass 
spectrum in Fig. [2j The BW width of each resonance is assumed to be constant and an 
overall three-body phase-space factor is applied. In the fit, the BW shapes are modified 
by the effective luminosity [l3| and m 7r + 7r -^,( 2S ')-dependent efficiency, which increases with 
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FIG. 2: The ir + ir~ip(2S) invariant mass distribution for events that pass the i^{2S) selection. The 
open histogram is the data while the shaded histogram is the normalized tp(2S) sidebands. The 
curves show the best fit with two coherent resonances together with a background term and the 
contribution from each component. The interference between the two resonances is not shown. 
The two dashed curves at each peak show the two solutions (see text). 
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FIG. 3: (a) The M r 2 ec and (b) polar angle distributions of the -K +r K~il)(2S) system in the e + e~ CM 
frame for the i: + it~iI)(2S) events with m^+^-^^s) ^ [4.0,5.5] GeV/c 2 . The points with error bars 
are data, and histograms are from MC simulation. 



m 7r + 7r -^,(25) from 3% at 4.3 GeV/c 2 to 5% at 4.7 GeV/c 2 . The effects of mass resolution, 
which is determined from MC simulation to be 3 MeV/c 2 -6 MeV/c 2 over the full mass range, 
are small compared with the widths of the observed structures, and therefore are neglected. 

Figure [2] shows the fit results with two solutions with equally good fit quality. In these two 
solutions, the masses and widths of the resonant structures are the same, but their partial 
widths to e + e~ and the relative phase between the two resonant structures are different (see 
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FIG. 4: 7r + 7r invariant mass distributions of events in different ■k + tt ip(2S) mass regions, (a): 
4.0 GeV/c 2 < m^+ v -^ {2S) < 4.5 GeV/c 2 , and (b): 4.5 GeV/c 2 < m^-^S) < 4.9 GeV/c 2 . 
Points with error bars are data while the histograms are MC simulation with the phase-space 
distribution generated at y/s = 4.4 GeV (a) and 4.7 GeV (b). 

TABLE I: Results of the fits to the Tr + ir~ifj(2S) invariant mass spectrum. The first errors are 
statistical and the second systematic. M, T tot , and B ■ T e + e - are the mass (in MeV/c 2 ), total 
width (in MeV/c 2 ), product of the branching fraction to 7r + ir~ip(2S) and the e + e~ partial width 
(in eV/c 2 ), respectively. <fi is the relative phase between the two resonances (in degrees). 



Parameters 



Solution I 



Solution II 



M(Y(4360)) 
r tot (Y(4360)) 



4361 ± 9 ± 9 
74 ± 15 ± 10 



S-r e+e -(y(4360)) 10.4 ±1.7 ±1.5 11.8 ±1.8 ±1.4 
M(Y (4660)) 4664 ± 11 ± 5 

r tot (y(4660)) 48 ±15 ±3 

B-r e+e -(y(4660)) 3.0 ±0.9 ±0.3 7.6 ±1.8 ±0.8 
6 39 ±30 ±22 -79 ±17 ±20 



Table HJ) |14j . The interference is constructive for one solution and destructive for the other. 
To determine the goodness of fit, we bin the data so that the expected number of events 
in a bin is at least seven and then calculate a x^/ndf = 4.7/3 corresponding to a C.L. of 
19%. The background level from the fit is 0.19 ± 0.14 events per 25 MeV/c 2 bin, in good 
agreement with the ip{2S) mass sideband estimate of 0.12 ± 0.05. The significance of each 
resonance is estimated by comparing the likelihood of fits with and without that resonance 
included. We obtain a statistical significance of more than 8a for the first peak (hereafter 
referred to as the K(4360)), and 5.8a for the second one (the F(4660)). 

The systematic errors in the mass and width measurements are dominated by the choice 
of parameterization of the resonances, especially the mass dependence of the widths; the 
range of changes in the fitted values for different parameterizations is reflected in the errors 
listed in Table [H Other sources of systematic error, such as the mass resolution and the 
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mass scale, are negligible. 

The uncertainties in B ■ T e + e - due to the choice of parameterization are 7% for the 
F(4360) and 10% or 3% for the two Y (4660) solutions. There are other sources of systematic 
errors for the B ■ T e + e - measurement. The particle ID uncertainty, measured using the 
e + e~ — > ip(2S) — > 7T + 7T~J/if) samples j3[, is 5.0%; the uncertainty in the tracking efficiency 
is 1%/track; the uncertainties in the J /if) mass, if>{2S) mass, and M r 2 ec requirements are also 
measured with a control sample of e + e~ — > if>(2S) — > 7i + 7i^J/if) events. For these events, 
the MC efficiency is found to be higher than in data by (4.3 ± 0.5)%; a correction factor is 
applied to the final results and 0.5% is taken as the associated systematic error. 

Belle measures luminosity with 1.4% precision while the uncertainty of the radiator in the 
PHOKHARA program is 0.1% 13j. The main remaining uncertainty in PHOKHARA 



is associated with the modelling of the 7T + ir~ mass spectrum. A MC simulation with 7t + tt~ 
invariant mass distributions that reflect the observations shown in Fig. H] yields an efficiency 
that is higher than the phase-space simulation by about 9%, which is used in the fits with 
half of the correction (4.5%) taken as the systematic error. According to the MC simulation, 
the trigger efficiency for the events surviving the selection criteria is around 98% with an 
uncertainty smaller than 1%. The uncertainty in the world average [l2| values for B(if)(2S) — > 
7T+7T- J/if)) is 1.9% and that of B{J/ip -> £+£-) = B{J/ip -> e+e") + B{J/if) -> is 1% 

where we have added the errors of e + e~ and fi + fi~ modes linearly. Finally, the statistical 
error in the efficiency is 1.3%. Treating each source as independent and adding them in 
quadrature, we obtain total systematic errors on B ■ T e + e - in the range 10-14% for the two 
solutions for the F(4360) and F(4660), see Table HH 



TABLE II: Systematic errors in the B ■ T e + e - measurement. 



Source 


Relative error (%) 


Parameterization 


3-10 


Particle ID 


5.0 


Tracking efficiency 


6 


J /if; mass, i/)(2S) mass, and M r 2 ec 


0.5 


Integrated luminosity 


1.4 


m 7T + n - distribution 


4.5 


Trigger efficiency 


1 


Branching fractions 


2.1 


MC statistics 


1.3 


Sum in quadrature 


10-14 



The cross section for e + e — > tt + 7t if)(2S) for each 7t + tt if)(2S) mass bin is calculated 
according to 

^obs _ n bkg 



e i C l B(if}(2S) -> n+iT-J/if))B(J/iP -> £+£-) ' 

where n° bs , £j, and Ci are the number of events observed in data, the efficiency, and the 
effective luminosity in the i-th 7r + 7r~if)(2S) mass bin, respectively; n bkg is the number of 
background events measured in ip{2S) sidebands, taken as 0.23 ±0.09 events per 50 MeV/c 2 
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for all the bins 
taken from Ref. 



B(ip(2S) -> tt+tt-J/^) = 31.8% and B(J/tp -> = 11.87% are 

The resulting cross sections are shown in Fig. [5], where the error bars 



16 . The 



include the statistical uncertainties in the signal and the background subtraction 
large error bars at low mass are due to the low efficiencies. The systematic error for the 
cross section measurement, which includes all the sources listed in Table HT1 except for that 
from the BW parameterization, is 9.5% and common to all the data points. 




cm 



(GeV) 



FIG. 5: The measured e + e — > 7r + 7r ip(2S) cross section for y/s = 4.0 GeV to 5.5 GeV. The errors 
are statistical only. Bins without entries have a central value of zero. 

In summary, the e + e~ — > 7T + 7r~ip(2S) cross section is measured from threshold up to 
5.5 GeV. The measured cross sections are consistent with results from BaBar Two 
distinct resonant structures are observed, one at m — 4361 ± 9 ± 9 MeV/c 2 with a width 
of 74 ± 15 ± 10 MeV/c 2 , consistent with the structure observed by BaBar in mass but 
with a much narrower width, another at m = 4664 ± 11 ± 5 MeV/c 2 with a width of 
48 ±15 ±3 MeV/c 2 , that has not been previously observed. The resonant structures reported 
here are distinct from the ones observed in e + e" — > 7r + 7r~JM 111, l3|j . There are no known 
vector charmonium states that match these measurements 12|, |l7|; according to potential 
model calculations 18, the 4 3 1 S 1 , 5 3 Si, and 3 3 -Di charmonium states are expected to 



be in the mass range close to the two resonances measured. We note that coupled- channel 
effects and rescattering of pairs of charmed mesons (£>(*)#(*), £>(*)£>(*)) may affect the above 
interpretation 
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